Quantum scalar fields in the half-line. 
A heat kernel/zeta function approach. 

J. Mateos Guilarte 1 , J. M. Murioz Castarieda 2 , and M. J. Senosiain 3 

1 Dcpartamento de Fisica Fundamental and IUFFyM, 
s * ■ Univcrsidad dc Salamanca 

o ■ 

f^*) ' 2 Dcpartamento de Fisica Tcdrica 

Universidad de Zaragoza, SPAIN. 
3 Dcpartamento de Matematicas 
\Q • Universidad de Salamanca, SPAIN. 



Oh- 



Abstract 



i -^h ■ In this paper we shall study vacuum fluctuations of a single scalar field with Dirichlet 

boundary conditions in a finite but very long line. The spectral heat kernel, the heat 
partition function and the spectral zeta function are calculated in terms of Riemann 
Theta functions, the error function, and hypergeometric pFq functions. 
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1 Introduction 

00 

CN I In collaboration with J. Sesma, J. Abad devoted part of the last years of his fertile scientific 
C^; career to studying the role of special functions in quantum field theory. In this brief memoir, 
elaborated to honor Julio's memory, we explore the influence of using Dirichlet boundary 
conditions in quantum field theory. Specifically, we shall address the Higgs model in (1 + 1)- 
dimensions but we shall restrict the spatial line to become a finite interval. Then, Dirichlet 
boundary conditions at the endpoints of the interval will be imposed on the field. Eventually, 



5_i ■ we shall allow the length of the interval to tend to infinity to describe the situation in which 



the mesons meet an impenetrable wall. Our playground is thus the analysis of scalar quantum 
fields living in a half-line. 

In this short work we shall concentrate on computing very basic quantities. Essentially, 
we shall deal with vacuum fluctuations in such a way that the spectral zeta function of the 
second-order differential operator governing small fluctuations around the vacuum will be used 
to regularize the divergent zero-point energy. The spectral information is also encoded in the 
associated A-heat partition function and A-heat kernel. These spectral functions permit a 
high-temperature asymptotic expansion, which, in turn, determines via the Mellin transform 
the meromorphic structure of the spectral zeta function in terms of the heat coefficients. The 
main sources of our approach are References [T], [2J, and [3] as well as [U] and [7j. We hope 
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that Julio would have been pleased with our results. In recent times he was one of those rare 
theorists trusted and praised by experimental and applied physicists. 



2 The Higgs model in a line 

In the (1 + l)-dimensional toy Higgs model the action 

S= [dy>{l^-^(y ,y)-^T 



2dy»d Vil 4 vr ^ u '* y A 

governs the dynamics of the scalar field ip(yo, y) : M 1 ' 1 — ► K. We choose the metric g^ v = 
diag(l, —1) in (l+l)-dimensional R 1 ' 1 Minkowskian space-time. In the natural system of units 
H — c — 1 the dimension of the field, the mass, and the coupling constant are respectively: 
[■0] = 1, [A] = [m 2 ] = L~ 2 . In terms of non-dimensional space-time coordinates and fields 

m VA 
the action functional and the field equations of the A (0)2 model read: 

m2 f J j / 1 ^ d( t> „\ ^2 



S = — / dx 2 < -TT^TT^ - -(r(x ,x) - 1) 

A J {2dx»dx fi 2 VV 
<9 2 d 2 (j) 

^|(x ,x) - — (x ,x) = 20(x o ,x)(l - 2 (x o ,x)) 

The shift of the scalar field from the homogeneous stable solution, <p{x^) = 1 + H(x tl ), leads 
to the action 
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m 



S = — l 



—dpHd^H — 2H 2 (x^) 



2H 3 (x^) + -H A (x^ 

which shows the spontaneous symmetry breakdown of the internal parity Z 2 symmetry. 



3 Zero point vacuum energy with Dirichlet boundary 
conditions 

The linearized field equations 

d 2 5H d 2 5H 

(x ,x) q^-(x ,x) +4SH(x ,x) = (1) 

allow us to expand the Higgs field 5H(x , x) as a linear superposition of solutions obtained by 
means of separation of variables: 

5H(x , x) = 1 • Hk)e-^°f(x; k) + a*(k)e ik °*°f (x; k)} . (2) 



([2]) is the general solution of ([T]) if the dispersion relation between the frequency and energy 
of the plane waves k$ — k 2 — 4 = (k — u>(k) = vF+4) holds. Of course, f(x; k) are the 
eigenfunctions of the second-order fluctuation operator: 

K = -^ + A , K f(x;k)=cu 2 (k)f(x;k) . (3) 

In the normalization interval I — [0,1], I = &t=, the spectrum of Kq with Dirichlet boundary 
conditions (following the method developed in [8]) 

K f n (x) = u 2 J n (x) , / n (0) = f n (l) = 

is: 

i 71 2 71-2 2 , a J- t \ . /7T s m+ 

k n = -n , u n = —n + 4 , f n (x) = W y sm(ynx) , rt G Z 

Therefore, the classical Hamiltonian is tantamount to an infinite number of oscillators given 
by the Fourier coefficients of these standing waves: 

rr(2) ™ 3 f j (IdSH d5H 185H 85H _ TT/ , rrT . 



oo 

m 



—j= io(k n ) ^a*(k n )a(k n ) + a(k n )a*(k ri 

* 71=1 



Canonical quantization [a(fc n ), a' (A; m )] = <5 nm promotes the Fourier coefficients to creation and 
annihilation operators and gives the free quantum Hamiltonian: 



^o 2) = ^f>(^) (a\k n )a{k n ) + \ 
It is clear that the vacuum (a){k n ) |0 >= 0,Vn) energy is not zero but: 



oo 

AE =< 0\H\0 >= ^=^(fc n ) = ^f Tr ^o 

* 71=1 * 

a divergent quantity. 

3.1 The heat function 

Better expectations of convergence are offered by another spectral function, the i^o-heat func- 
tion: ^ ^ 

Tr D e-P Ko =[ dxK Ko (x,x;(3) = f]e~' 3 ^ +4 ') (4) 
where K Ko (x, y; (3) is the kernel of the i^-heat equation 

A + R^j m{(3,x)=(^-^+ 4^ Kxo (z, y; P) = , tf* (x, y; 0) = 5(x - y) 



and (3 = is proportional to the inverse temperature. Moreover, via the Mellin transform 
the spectral zeta function is obtained: 

( K Js) = — - / dpp^Trne-^ = V — ^ . (5) 

^ KoK } T(s)J Z^^ + A y ^ 

We shall use this meromorphic function of the complex variable s (and will return to this later) 
to regularize the divergent sum of vacuum fluctuations, AE , by assigning to it the value of 
the series at a regular point in the s complex plane. 



3.1.1 Riemann Theta constants 

The .Ko-heat function is essentially given by a Riemann Theta constant: 



TrDe -^o = ^ e -/ 



2 e 



-4/3 



J2 exp[-/3^n 2 ]-i 



n=l 



-4/3 












7T 



(Ojz — /3) - 1 



Here, we denote the very well known Riemann or Jacobi Theta functions in the form: 



(6) 







a 
b 



(z\r) = ex P 



2m[(n + a)(z + b) + -(n + afr] 



z G C , r G C, Imr > , a,b — 0, \. Thus, we need the Riemann Theta function at the 
2 = point (Theta constant), the modular parameter r = ip [3 (determined by j3 and /), and 
the "characteristics" a = b = 0. Use of the Poisson formula 












(0,4/3) = 4^0 








(0,i 



7T/3' 



allows us to write the ii'o-heat function in the new form: 



-4/3 / l 












(0,i 



7T/3' 



From this, an asymptotic formula for the behavior of the i^o-heat function is obtained: 












(0,z-^)-l+O(e-f; 



-4/3 / l 



/3^0 



1 )+0(e 



(7) 



3.1.2 Physicists' derivation: the Error function 

We now offer a derivation of the asymptotic formula by means of physicists' techniques. The 
idea is to look at the problem when / is very large: I — > 00. The spectral density of the 
standing waves can be determined from the phase shifts 5 D {k) = —Si(2kl) (Si(x) is the sine 
integral function) due to the reflected waves: 



sin(kl + 5 D (k)) = 



kl + 5 u {k) =mr , fi£Z 



1 



*w = ^ = -+e-^) = -(i 



dn I 1 dS 



D 



sin 



(2kl) \ 



kl J 



dk 71 7T dk ' 7T 

Thus, we end with an integral, rather than a series, for the i^o-heat function in terms of the 
error function: 



T -?k = L , dk i x 



7T 



A'/ 



7 



Erf 



VP. 



(8) 



The high-temperature formula agrees perfectly with (J7J) 



Erf[-^=] =^ l + C»(e-f; 



TV rvj 

lr_D e =/3^o 



-4/3 / j 



C(e" 



and, neglecting exponentially small contributions, we find the coefficients of the high-temperature 
expansion: 



-4(3 



' ±Erf[ ' 



4tt 2 



z 



c (K 



, Ci/ 2 (K ) 



Cn(K ) = 



Vra > 1 . 



3.2 The spectral zeta function 

3.2.1 Epstein zeta function 

Mellin's transform of the i^o-heat function ([6]) provides the spectral zeta function in terms of 
the Epstein zeta function E(s,a\A) = 



i 



n=— oo ( J 4n 2 +a) s ' 



1 T 00 1 Z" 00 / °° 2 



1 °° 1 

i V - 

9 / -> f^ 2 2 



1 1 , ,7T 2 , 



2 ^ (nl n 2 + 4: )s 2 2s + 1 2 v ' 'Z 2 ' 2 2s +! 

n=— oo V ' 

Mellin's transform, however, of the Poisson inverted version 



V W WV nSZ/{0} 

gives the spectral zeta function as a series of modified Bessel functions of the second type. 
Moreover, formula (Q shows that there are poles of (k ( s ) a ^ the points 

= I_I_3_5_Z 2j + l + 

s 2 , 2 , 2 , 2 , 2 ,---, 2 Jfc^ 

because -Ki/2_ s (4Zn) are transcendental entire functions, i.e. holomorhic functions of s in C/oo 
with an essential singularity at s = oo. 
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3.2.2 Physicists' derivation: Hypergeometric p Fq functions 

Mellin's transform of the (jHj) version of the i^o-heat function 

I 



_L r( s -i/2)/- 1 



• v 



F 2 [s;l/2 + s,l + s;4l 2 }) , (10) 



supplies a third analytical expression of the spectral zeta function. Euler T functions and 
hypergeometric pFq functions, with power expansion around z = 



pFq[cii, a 2 , ■ ■ ■ , a p ; 6i, 62, • • • , z] = 7^ 



:0 (&i)fc(& 2 )k • • ■ (b q ) k ' fc! ' 

where (a)k = a(a+l)(a + 2) ■ • ■ (a-\-k — 1) is the Pochhammer symbol, enter the third formula 
of Ckq( s )- It is clear that the physical point s = — | is a pole of at least T(s — |). Other 
poles come from the other poles of V(s — \), s — | = 0, —2, —3, —4, —5, • • •, and the poles of 
iF 2 [l/2; 3/2, 3/2 — s; 4/ 2 ] and i-F 2 [s; 1/2 + s, 1 + s; 4/ 2 ], which are meromorphic functions of s. 
From the residue representation of these functions 

pr i 3 3 r(pr(f-s) ^ / r(i-^)(-4/ 2 r r , A , 

\to -; -, s; 4/ = — - — ■ > res u — - — : — jr. -1 (u) (—7) 

2 2' 2 ' J r(|) f^ Vr(| - u)r(| - s - u) y - ) ) yj > 

w r 1 , , 4/2l r(| + 8 )r(i + 3 ) ~ / r( g - M )(-4/ 2 )- , 

iF 2 s; - + s, 1 + s; Al\ = — > res u —j— — — — -T{u / 

2 r(s) V r (i + s - u)T{l + s - u) 

we find poles when 3/2 — s = —k±, 1/2 + s = — fc 2 , 1 + s = — k 3 , fci, fc 2 , &3 G Z + |J{0}- All 
together, there are poles of Ck ( s ) 

s = . . • - 5/2, -2, -3/2, -1, -1/2, 1/2, 3/2, 5/2, 7/2, 
3.3 The heat equation kernel 

Finally, in this sub-Section we analyze how the i^o-heat function, henceforth the spectral zeta 
function, are obtained from the i^o-heat kernel. 

3.3.1 Jacobi Theta functions 

The .Ko-heat equation kernel satisfying the Dirichlet boundary conditions 

+ 4)Kj> (x,y;P) = , Kj> o (x,y;0) = 5(x-y) , AT£ o (0, y; (3) = K° Q (l, y; (3) = 

(11) 



dp dx 2 
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is: 



K Ko( x >V>P) = ye 4/3 J^ sin(ynx)sin(yn?/) • e 



n=l 
-4/3 00 



zd„2 e 



-4/3 



S y^ j sin(— nx)sin(— ny) • e 



2/ 



^cos(y?i(x — y)) — cos(yn(x + y)) j • e /3 ^ T ™ 



-4/3 



2/ 

















(12) 



Alternatively, a modular transformation allows us to express the heat kernel in the new form: 



-4/3 



■ e 4 f 6 








.j X -y . . I 2 , 

2/3 'V/T 



_ e 4/3 








.,1 + y . . / 2 
-iL z - U- 



2/5 1 7T/3' 

because the Jacobi theta functions involved are modular forms of weight 1/2. 

3.3.2 Physicists' derivation: the Laplace transform 

Another route to solve pip is to look for solutions of the form 

K° Q (x, y- (3) = K Ko (x, y;(3) + e" 4/3 J D(x, y; /?) 

where 



(13) 



K Ko (x,y;f3) 



-4/3 



cxp [ 



(x - y)' 



V^P 4/5 1 

is the .Ko-heat equation kernel with periodic boundary conditions. (TT3]) complies with Dirichlet 
boundary conditions if: 



D{x,y;0) = , D(0,y;P) 



1 



e~V , D{l,y;/3) 



1 



e ^ . (14) 



The Dirichlet boundary conditions ( Tl4l) force the Laplace transform of D(x, y\ (3), D(x, y; s) 
J df3 e~ s(3 D(x, y; (3), to satisfy: 



D(0,y;s) 



D(l,y; s) 



-sT<i-y) 



Moreover, the ansatz (TIB"]) solves (TTTT) if D(x,y; s) solves the Laplace equation: 



(15) 



(16) 



The general solution of (fT6l) is 



D(s,j/; S ) = A( 2 /)e-^ a! + 5(j/) e V Sl 
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which complies with (fT5|) if: 



D(x,y; s) 



2^~s 



+ 



exp[-y/s(Z + x - y)] - exp[-y/s(x + y - I)] 

ply/s g — lyfs 

exp[-A/i(/ - x + y)] - expj-y^ - x - y)} 

ply/S p — l^JS 



;i7) 



The last step is to take the inverse Laplace transform of D(x, y; s) as given in (j!7p . To do 
this, it is convenient to write the common denominator as a power series expansion: 



1 



-ly/s 



E^ 



-(2n+l)ly/s 



n=0 



or, 



D(x, y; s) = — = ■ ^ [exp[-v^(2i(n + l) + x-y)]- exp[-Vs(2nZ + x + y)] 

* n=0 

+ exp[-v^(2Z(n + 1) - x + y)] - exp[-^(2l(n + 1) - x - y)]] 
The inverse Laplace transform of this is easy and gives: 
1 



D{x,y-(3) 



V^P 



n=0 



. (2l(n + 1) + x - yf (2ln + x + y) 2 1 
exp[ — J — exp[ 



4/3 



4/3 



r (2l(n + l)-x + y)\ {2l{n + l)-x-yf 
+ ex P[ T7, J-exp[ 



4/3 1 4/3 

From this formula we derive the Dirichlet Ko-heat kernel at coinciding points 

-4/3 



K° (x,x;P) 



n=0 



o _ e /3 - e 



which in turn provide the i^o heat function through integration on the interval: 

Tr D e-^° = 



l /p-4/3 

dxKl(x,x;P) = ^==[l + 2j2(-lT 



V^ 



e 



n=0 



-4(3 



V^P 



v^p 



n=0 



2-0 





1/2 



dx { e v + e 



E(- X ) n ( Erf 



n=0 



V ' 7T/3 ; 

Z(n+1) 
VP 



Erf 



Because 



£(-!)" (Erf 



n=0 



Z(w + 1) 



Erf 



In 

VP. 



.VP\ 



we again find 

in the high-temperature regime. 

4 Summary and outlook 

In sum, we have found three different expressions for the i^o-heat function: 

p-4/3 _ p-4/9 -4/3 



where 




Figure 1: Plot of: a) /xflrl), b)/ 3 (|r|), and c) / 3 (|r|). 

Figure 1 shows the Mathematica graphics of /i(|r|), /2(|t|) and /3(|r|). In Figure 2(a) 
the graphics of /i(|t|) and /2(|t|) are shown together. Simili modo, the graphics of fi(\r\ and 
/3(|t|) are plotted together in Figure 2(b). It is clear that all three graphics agree perfectly 
when (3 — > (high-temperature) and/or I — > oo (infinite length of the interval). /i(|t|) and 
/2(|t|), however, start to differ at |r| = 0.7, whereas there are no differences in the graphics of 
/i(|t|) and /s(|t|). It is amazing how two different derivations involving highly sophisticated 
special functions lead to identical curves ! From a physical point of view we are tempted to 
speculate that /3(|r|) would give the exact result /i(|t|) because the infinite rebounds of the 
standing waves in the walls at x = and x = I are accounted for. Instead, /^(a?) counts a 
single rebound in the x = wall, which is a legitimate approximation for I — > oo. 
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f_l,f_2 f_l , f_3 




Figure 2: Plot of: (left) /i(|t|) (continuous line) and /2(|t|) (dashed line), and (right) /i(|r|) 
(continuous line) and /3(|r|) (dashed line). 

We plan to follow this work by extending these computations to the kink sector of the 
model. The idea is to compute the one-loop kink mass shift in the framework developed in 
Reference [5] using Dirichlet boundary conditions instead of the periodic boundary conditions 
that are more conventional in quantum field theory . It will also be of great interest to perform 
the same program using more general families of boundary conditions, combining the method 
developed in [H [5] with the formalism developed in references [HI EH IE]- 
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